Contract  SD“88 


Tochnicit  rspert  Nc  A 8  PA  "15 

PERTURDATION  THEORY  OF 
THE  SELF-CONSISTENT  FIELO 


A.  1C.  fiajagapiJ 


Cl'cARlNGHOUSE 
FEI.'ERAI>  SC^ENriFIC  AND 
TECHmCAL  INTORMA-RIOK 
Mltrefiehel  * 


9 

L-X 


A  orA  %O.Sd  UA  ppL%| 


L 


mm 


DDC 


May  20, 1SSS 


DOCiRA  F. 


OmSIOfl!  OF  Ef^GJ^iEERINC 
HARVARD  U’^IVFRSITY  • 


mm  PHYSIC 
.  MASSACHi 


PERTURBATION  THEORY  OF  THE  SELF-CONSISTENT  FIELD 


by 

A.  K.  Rajagopal 


Technical  Report  No.  ARPA-15 
Contract  SD-88 
May  26,  1965 


Submitted  to: 

Advanced  Research  Piojects  Agency 
The  Department  of  Defense 


tion  of  this  document  is  unlimited. 


Division  of  Engineering  and  Applied  Physics 
Harvard  University 
Cambridge,  Massachusetts 


ARPA- 1 5 


ABSTRACT 


A  perturbation  theory  of  the  eelf -conaistent  field  is  here  developed 
using  a  Green's  function  method  after  the  manner  of  Peng.  It  is  shown 
that  a  self-consistent  perturbation  theory  yields  the  same  answers  as  the 
conventional  random  phase  approximation.  The  equivalence  of  the  perturbation 
approach  to  spin  waves  in  metals  developed  by  Herring  and  the  random  phase 
approximation  for  the  same  problem  is  seen  to  be  an  example  of  this.  The 
effect  of  the  weak  periodic  field  on  a  Hartree-Fock  electron  system  is  here 
explicitly  done  both  by  the  perturbation  approach  and  the  the  random  phase 
approximation  with  the  Green's  function  approach.  The  Green's  function 
approach  is  finally  generalized  so  as  to  include,  in  principle,  correlations 
beyond  the  Hartree-Fock. 


PERTURBATION  THEORY  OF  THE  SELF-CONSISTENT  FIELD 


by 

♦ 

A.  K.  Rajagopal 

Division  of  Engineering  and  Applied  Physics 
Harvard  University,  Cambridge,  Massachusetts 


INTRODUCTION 

About  twenty  years  ago,  Peng  wrote  a  paper  [1]  with  the  above  title, 
developing  a  straightforward  perturbation  theory  of  the  Hartree-Fock  (HF) 
self-consistent  field  equations  when  there  is  a  small  perturbing  term  in  the 
Hamiltonian.  This  method  was  used  a  decade  later  by  Herring  [2]  in 
demonstrating  the  existence  of  spin  waves  in  metals.  Recently,  this  problem 
has  been  reexamined  by  many  authors  in  the  random  phase  approximations  (RPA). 
In  particular,  the  present  author  noted  [3]  (where  references  to  other  work  could 
be  foimd)  the  close  relationship  between  Herring's  work  and  the  RPA  results  he 
obtained  by  a  Green's  function  approach.  But  at  the  time  of  writing  about  it, 
we  were  not  aware  of  an  explicit  relationship  between  the  two  apparently 
different  approaches  to  such  problems,  although  the  equivalence  of  the 
two  methods  was  surmised.  In  this  report,  a  proof  of  this  equivalence  is 
explicitly  constructed,  in  response  to  repeated  demands  for  one  by  Professor 
Brooks  earlier  and  more  recently  by  Dr.  Herring.  It  might  also  be  of 
interest  to  mention  here  that  the  work  of  Thouless  [4]  on  RPA  is  closely 
related  to  that  of  Peng  [1],  even  though  Thouless  uses  a  variational  method 
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and  Peng,  a  perturbation  approach.  Slater  (5],  in  his  Quantum  Theory  of 
Atomic  Structure  .  Voliime  I,  has  a  demonstration  of  the  following  simple 
theorem:  The  equations  obtained  in  a  variational  principle  involving  a  certain 
basis  ef.t  of  functions  4re  the  same  as  those  involving  a  perturbation  approach 
where  the  same  basis  set  is  used  as  the  zero-order  functions.  The  present 
report  demonstrates  this  even  when  one  has  nonlocal  self -consistent  field 
e  ':)uations . 

In  most  many-electron  calculations,  one  deals  with  a  uniform  electron 
gas  with  a  uniform  compensating  positive  background.  We  first  treat,  by  the 
present  techniques,  the  effect  of  a  weak  periodic  positive  background  on  a  HF 
electron  gas.  (Professor  Brooks  informed  the  author  that  he  gave  this  problem 
in  one  of  his  advanced  courses  on  quantiun  theory  of  solids  ten  years  ago  at 
Harvard  University.' )  This  is  essentially  the  nearly  free-electron  problem 
when  we  have  HF  correlations.  See  Anderson  [6]  for  a  discussion  of  the  same 
in  the  ab8en-:e  of  HF  correlations.  We  do  this  problem  first  by  the  Peng 
method  and  then  by  the  Green's  function  approach.  Finally,  we  generalize 
this  Green's  function  approach  formally  when  the  electron  interactions  are 
taken  into  account  beyond  the  HF  approximation.  This  problem  was  suggested 
to  the  author  by  Professor  Ehrenreich.  The  results  are  quite  similar  to 
the  recent  work  df  Tanaka  [7],  who  uses  a  diagrammatic  approach  to  this 
problem. 
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PENG'S  PERTURBATION  EQUATIONS  AND 
RANDOM  PHASE  APPROXIMATION 
We  first  collect  Peng's  [1]  results  a.nd  compare  them  with  the  RPA  of 
Thouless  [4],  Then  we  go  on  to  develop  &  Green's  function  reformulation  of  the 
same  type  of  perturbation  theory  in  the  HF  approximation.  We  then  explicitly 
demonstrate  the  equivalence  of  our  equations  obtained  here  with  those  obtained 
earlier  [3]  for  the  case  of  spin  waves  in  an  electron  gas.  In  the  next  section, 
we  examine  the  effect  of  a  w'eak  periodic  field  on  a  many-electron  system. 

Peng  proceeds  as  follows.  The  HF  equations  for  an  n>c'ectron  system  are 
(Tt  V  +  G  -  E^)  =  0  (T  =  1,  2,  . .  .  n)  (1) 


T  is  the  one-electron  operator  containing  the  kinetic  and  potential  energies. 

V  is  the  one-electron  external  potential  to  be  treated  in  perturbation  theory. 
(Peng  used  a  moregenexal  set  of  perturbing  terms;  we  have  simplified  it  here 
for  the  present  purposes. )  And  G  is  the  HF  operator  defined  by 
n 

\=1 

the  first  being  the  Coulomb  term,  and  the  second  the  exchange  term,  and 


(3) 


where  v  is  the  two-particle  interaction  potential  which  is  self-adjoint  and 
symmetric.  s  is  the  spin  variable.  In  most  cases,  v  =  e^/|  r  -  r'  { .  The 
energy  parameters  E^  are  the  HF  one-electron  energies, 
unperturbed  HF  problem  be  assumed  to  have  been  solved: 


Let  the 
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(T  +  G  -  ^  0 


(4) 


where 


(5) 


X.*l 


Note  that  iiince  the  operator  G  is  self  adjoint  and  line,  r,  the  set 
is  orthoaorTnal.  One  may  therefore  write 


py 


Writing 


ap 


P 


ap 


and 


(A  (®>  (A  «»  V  A  <«> )  A  <«> )  =  V 

'^a  ’  ^p  '  aftvp 

Eq.  (1)  may  be  written  as 


J 


(6) 


(7) 


I 

P 


T  +  V  + 

Qp  ap 


nr 

\=1  pv 

(1) 


li\  ^'^Q/Avp  '^QjLipv' ’^vX,  ■“7''ap|  ^P7 


^R'v  ~  ^ 


Let  us  denote  V^p  =  V^p  ,  being  of  first  order  in  perturbation. 
Then  write 

^P7  ^p-y  ^P7  ^  By  •• 


(8) 

(9) 


Py 


E  =  E  E  E  ^^^  +  .  .  .  J 

y  y  y  y 


(10) 
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One  then  has 


and 


n 


E  T  +  y  (v  . .  -  V  ) 

a  aa  ^  aX.Xa  aXaX 


X=1 


n 


ay  oXKy  a\y\  \  r 


X.=  l 


J 


(11a) 


and 


and 


n 


X=l  n 


^aX^a  flX  *  ^'^Q/iXa~ ^ 


(11b) 


(E  E  )*?/ y  y [(v  -V  -V  )‘2/ ] 

a  y  ^ay  /_,  akfiy  aKyii’^nX.  ofiky  ^aLtiy\  ^fik  ■* 


X=1  ii 


'a  "4 


X=1  (1 


*liXa  allay ‘^nX  J  + 


\ 


n 


+y  y(v  -V  )^ 

L  L  ativa  ofMav’^nX  ^  v\ 


X=1  liv 


n 


X=1  ^ 


V.  -V,  )Y["’  + 

a\tiy  akyfi  ~  tik 


^^a4iky~  ^anyk'  nk 


)  P>*  ]  = 


n 


=  -  (E  ^^^-E  y  (v  -v  )2<^  ia^y) 

a  y  ^  ay  LiLi  a^tvy  atiyv  ^tik  vk  '  r  t/ 


X=1  nv 


J 

(He) 
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Moreover,  the  normalization  condition  gives 


2(^(2)  iz)*  ^  \  (1) 

p 


(lid) 


etc 


Also  nuite  generally  one  has 

♦ 

^Pvna  ~  '^ativP  ~ 


(lie) 


The  homogeneous  counterpart  of  the  second  of  the  equations  (11b)  will  now 
be  identified  with  the  time>independent  RPA  equations  of  Thouless  [4]  . 

Svq>pose  that  creates  a  particle  with  wave  function  i^.  .  The  state 

v^ith  n> electrons  occupied, is  the  Slater  determinant 


I  «o  ^  °  <  11  '^i'^  •  I  0  > 


(12) 


i  =  l 

with  |0  >  representing  the  vacuum  state  and  such  that  ^ 


Thouless  [4]  showed  that  any  n- particle  Slater  determinant  |  which  is 


not  orthogonal  to  |  can  be  written  i.i  the  form 


!  »>  = 


exp  ^  ^  (f  ^ 

^  mi“m 

L  i=l  m=n+l 


«  > 

o 


(13) 


where  are  uniquely  determined  and  conversely  that  any  wave  function  of 


this  form  is  an  n- particle  Slater  determinant.  Now  suppose  I  ^  1®  fbe  HF 
state.  Choose  so  as  to  minimize  the  average  energy:  E  =  < #| H | #>/< ^| ®>. 


This  gives  eqviations  of  the  form  (11b)  satisfied  by  ^ \  a-nd  these  are  indeed 
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the  RPA  equations.  Note  that  to  all  r.x'ders  the  homogeneous  parts  of  the 
perturbation  eqiiations  are  the  same  and  of  the  form  of  the  homogeneous  part 
of  (11b).  When  we  generalize  this  in  the  Green’s  function  formalism  to  a  time- 
dependent  fonri,  we  see  the  identification  with  RPA  equations  even  more  closely. 

Both  the  Thouless  and  Peng  calculations  employ  as  starting  point  the 
complete  set  of  the  unperturbed  HF  solutions.  In  ordinarv  Qti^^ntxim  mechanical 
problems,  Slater  [5]  showed  that  in  such  a  procedure  the  perturbation  equations 
and  the  equations  obtained  from  the  variational  methods  are  identical.  This  is 
seen  to  be  the  case  even  for  the  equations  of  the  HF  type.  Recently,  Fukuda[8], 
without  the  knowledge  of  the  work  of  Thouless  [4],  showed  that  whenever  one 
has  r  ;rturbations  of  the  HF  states  due  to  any  external  fields,  the  new  state  can 
be  written  in  the  form  of  a  unitary  transformation  of  the  type  (12).  Based  on 
this,  he  constructs  stability  conditions  for  the  HF  state  when  one  has  v'arious 
external  perturbations.  This  again  led  him  to  RPA  equations,  again  in 
conformity  with  Peng's  calculations. 

We  will  now  develop  a  simple  Green's  function  formulation  of  the  same 
basic  ideas  as  in  Peng's  work  except  writing  it  in  a  time -dependent  form  which 
has  some  advantages.  L«t  us  consider  the  HF  Green's  function  equation  [3] 
in  the  presence  of  an  external  perturbing  field  V  (1). 


-U(l)- V(l)  +  i  j*d3  v(l-5)tr 


^HF^'^)- 


dl  v(l-3) 


Gj^p(ll)  =  6^^^(l-2) 


(14) 
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wh«re  G(12)  =  -p  <  T  (t// (1 )  (2))>  in  this  usual  notation  and  1  stands  for 

the  space-time  point  (r^tj).  Here  U(i)  is  the  one-particle  potential  in 
whos*  ;r':sence  the  problem  with  V  =  0  is  assumed  to  be  solvable.  We  have 
the  completeness  relation 


r  (i  2)  Gj^j,(2  1‘)d2  =  ^ 

> 

J  Gjjpdl)  G^jT^Il'ldl  6^^^(1-1') 

When  V  ~  0  ,  the  zero-order  HF  eqiuition  is 


dJv(l-3)  tr 


ChF  - 


-‘I 


d3  /(1.3)Gj^^°^(l  3)  Gj^^®^3  2)  -  6^^^(l-2) 


(15) 


(l^) 


and 

2)  G^p°>  *(Il')dr  = 


(17) 


We  now  develop  a  perturbation  theory  in  V.  For  this  purpose  we  write 


6G„^(12) 


v(T)dJt^jy 


6^Gj^p.{12) 


V(I,)V(l2)dTjdi; 


and 


.r 


aG^(i2) 

a  V  d) 


V(5)d5+  ... 


(18  ) 

{18’) 
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First  let  us  use  the  relation  (15)  to  write  £q.  (14)  in  the  form 

,2 


V  ^ 

i||“+  -  U(l)-V(l)+i  yd5v(l  -3)  tr  G^p(J3^ ) 


5<^)(1.2). 


.iv(l-2)Gj^p(12)  . 

Then,  collecting  equations  to  each  order,  we  have 


(19) 


*(12)= 


i  |^+  ^  -  U(  1 )  +  ijd Tv(  1  - J)  tr 


fi^^\l-2)-iv(l-2)Gjj^^l2) 

(20a) 


fiG  "p(12) 

—  ~  *  6(1-2)6;T-3)  f  i 


i  Jd5v(l-3)  tr -  6^^^U2)  - 


-  iv(U2) 


6G„p(12) 

. 


etc 


(20b) 


Now  from  (18),  (18'),  (15)  and  (17)  we  have 

OG  (12)  p  .  .  OG  '^23)  ... 

"  *  J  - 5^TTr~  Gj^'(^2)d2  dl  . 


T7W 


SV(3)  ^'HF 


(20c) 


■  1, 


dG^Q*(i2) 

Substituting  this  in  (20b), we  get  an  equation  for  — -  to  be  solved: 


6Gj^y(12) 
«  VT3)-' 


=  -  6^^^(l-2)6^^^(l-3)-i^ d2d3d'?v(l-3)t: 


(20d) 
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Thl*  ia  just  the  RPA  equation  for  [3] .  Note  the  appearance  of  Gj|p  in 
thia  equation.  It  ia  clear  that  we  can  generalize  the  approach  when  we  have  a 
more  general  Once  (20d}  ia  solved, one  puts  it  back  in  (18)  and  looks  for 

the  pole  in  the  energy  plane  in  G.  We  do  this  explicitly  for  the  case  of  the 
weak  periodic  field  in  the  next  section. 

At  this  point  we  would  like  to  draw  attention  to  the  equivalence  of 
Herring's  perturbation  approach  [2]  and  RPA  [3]  in  the  problem  of  spin  waves 
in  metals.  Herring  calculated,  using  Peng's  equations,  the  spin  wave  energy 
by  applying  an  external  static  transverse  field  to  second  order,  it  being  a 
^ime-independent  perturbation  calculation.  We  [3]  studied  the  same  problem 
by  an  approach  similar  to  the  one  outlined  above  in  the  last  paragraph  using 
a  time-dependent  exterxxal  field,  and  using  RPA.  The  eigenvalue  of  this 
equation  was  the  same  as  the  spin  wave  energy  obtained  by  Herring.  This 
answers  the  query  raised  by  Brooks  and  Herring  in  private  communication. 

In  the  next  section,  we  work  out  an  explicit  example  of  this. 


WEAK  PERIODIC  FIELD  IN  A  MANY-ELECTRON  SYSTEM 
We  first  study  this  problem  by  the  Peng  approach  and  then  derive  the 
same  resmts  and  generalizations  thereof  by  the  Green's  function  method. 

Let  the  weak  periodic  field  be  V  (  r )  such  that 


V 


K 


giK  r  reciprocal  lattice  vectors). 


(21) 
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In  the  absence  of  V(  r}f  the  Hi'  equations  of  a  uniform  electron  gas  in  a 
uniform  positive  background  admit  of  plane  waves  as  solutions.  So,  we 
employ  these  as  our  basis  when  V  occurs  in  the  problem.  Thus, 

y-  i(k+K,).r 

*  X(k  +  Kj)  .  (22) 

k" 

1 

Then  the  resulting  equation  for  X , proceeding  as  in  the  last  section,  is: 


X(k  +  K) 


(Kj)  X(k  +  K-Kj) 


(26) 
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The  first-order  equation  gives 

E^^\k)  =  V^^^(O)  =  0  for  simplicity 


(27a) 


and 


{E^®^(k+K)-E^®''(k))X^^\k+K)+^  [v(-K)-v(k-kj+2K)]  (  X^^^{kj+K)+ -K))  = 


“1 


=  -  V^^^(K) 

Neglecting  the  exchange  term  v(k-kj+2K),  we  find  that 


(27b) 


X^^^(k  +  K)  =  - 


V^^^(K)  1 


E^°'(k  +  K)-E^°^k) 


(27c) 


where 


€  (K)  =  RPA  dielectric  constant 

=  1  +  2v  (K)  ^ 

k  (occ) 


1 

E^^^(k  +  K)  -E^^\k) 


(27d) 


Note  the  appearance  of  HF  energies  everywhere  .  The  second-order  energy 
shift  is 


E^^^(k)  =y  V^^^(Kj)X^^\k-Kj)+y  y[v(-Kj)-v(k-kj+Kj)]  X^^^kj+Kj)X^^^k+Kj) 


K, 


l-j 

k,  K, 


+  y  y  [v(K,)-v(k-k,)l  X****(kj-Kj)  x“’(k^K^) 


k,  K, 


-  ' _ _ 
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which  simplifiea  to 


-I 


V*‘*(K  ) 

\-7WT~ 


E'  '(k+Kj)  -  E""(k) 


after  neglecting  the  exchange  terms  v(k-kj+Kj)  and  v(k-kj).  Note  that  the 
weak  periodic  field  gets  screened  by  the  RPA  dielectric  constant. 

Let  us  now  do  this  problem  using  the  Green's  function  approach  but 
including  formally  more  general  correlations  than  HF.  Specifically,  we 
rederive  {28)  ar  a  byproduct.  To  this  end,  we  have  to  solve  the  foliov.ng  set 
of  easily  derived  eqijations  which  are  all  exact: 

i|_+^  .  G(ll')  +  i  J  v( i-T)  G (ll) r(l4;l)K (IT)  GC-Tl'ldldldld? 

=  (29) 

"(12;  3)  =  -2)  a^^^{l-3)  -  i  Jv(l  -T)  G  (1  4)  F  (??;  3)  G  (SI)  Fd 2;  I)  K  (IT) 

dT  dl  d  J  d?  dl 

+  iJv(l.r)G(l  2)  -^^~^^K{ir)  dl  d2  dl  + 

+  iJv(l-T)G(l  2)  F(2  2;1)  dl  d2  dl  (30) 

K(12)  =  6^^^(l-2)  +  i  J  v(l-r)G(T2)  F(2l;4)  K(42)  G(3T'*^)  dT  d2  dl  d4  (31) 


(32) 


ARPA-15 


*14* 


where  is  introduced  only  to  generate  the  proper  terms  arising  from  the 

Interactions  between  the  particles.  contains  the  weak  periodic  field, 

^ext  ’  Coulomb  self  energy.  We  now  make  ofome  approximations. 

We  take 

r^®^(12,'3}  ^  .  (33) 

It  must  be  stressed  here  that  this  is  not  the  usual  HF  approximation.  Then 


.(G) 


(12)=  G(ll)  K^^'^dZ)  G{Il'*’)  dT  dl  (34) 


and 


“1 


V  ^ 

.  i  .  V 

|^8tj  2m  eff 


(1) 


v(1-1)G(12)K(2T)G(2  r)dl  dl  (35) 


=  6^^^(1-1') 


Going  to  the  Fourier  representation  we  get 


“-5^  C(kk').^V^jj(kj)G(k-kiik')- 
/  >^1 

-iT  vikj)  GSk+kpk^+kj)  K(k2;kj)  G(kjk')  =  6|^,  (36) 

kjk'^kj 


and 


Kskk')  =  -i  v(k)  2  G(k+k,;k2+k3)  Glkji  kj)  K(k3k')  .  (37) 

kik2k3 

In  the  absence  of  the  periodic  potential  and 
G'kk')  =  G^(k)0^, 


we  have 
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being  a  uniform  system,  so  that 

®  kk' 

K(kk')  = 

with 


I 


e  (k)  =  1  -  iv(k)  >  G^(k  +  kj)  G^(kj) 


(38) 


'1 


and  (36)  takes  the  form 


4  -  ‘2,  ^  G<«>(k  +  k,)  G(kk')  V^„(k,)  G  (k-k,;  k') 
k, 


Now 


1 


V^„(k)  =  J  K(k; 
K 


kk'  . 


K'  V^^^(K) 


(39) 


using  only  the  periodic  field,  so  that  one  finally  has 


G(kk')-^  -^|^^G(k-K,  k‘)  = 
K 


(40) 


where 


M(k)  =  mass  operator 


.T 

=  TTITT 

‘'i 


G^®^(k  +  kj) 


(41) 


This  is  the  dynamically  screened  exchange  interaction,  (k^  here  stands  for 

the  four  vector. )  Note  also  that  V^^^(K)  is  statistically  screened  by  c(K), 

(0) 

the  generalized  dielectric  constant  involving  the  G  '  which  is  not  the  HFG. 
Thus, 
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“  ‘  2m 


K, 


V^^^{K,) 

-^^G(k-Ki;k') 


(40') 


and 


u 


k-K 


TST 


“1  Y"  ^K,) 

-M(k-Kw)|  G(k-Kjk’):^)  ^  G(k-K-K^,k') 

~(1) 

4^  o,». 


(42) 


80  that  one  finally  has 


,u  - 


M(ku) 


■I 

K, 


TKTT 


(k-K,)^ 

u,-- ^-M(k-Kj;u,) 


G(kk‘) 


= 


kk' 


(43) 


The  po^es  of  G(kk')  determine  the  single  particle  behavior.  Writing  the 
unperturbed  sin^jie  particle  energy  in  the  form 

E<“>(k)  =  ^  +  M(ki  ) 


(44) 


the  modified  energy  dv  -o  weak  periodic  field  away  from  the  Erilluoin  zone 


boundaries  may  be  v.'ritten  as 


u  =  E^°\k)-  ^ 

K, 


1 

E^^^k-Kj)  -  E^°^(k) 


(45) 


If 


M(k;  ^  )  =  M^j,(k)  ,  €(K)  =l-riv(K)^Gj^^°^K+kj)Gj^^°^(kj) 
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we  recovei  the  residt  obtained  earlier  as  expression  (28).  At  the  zone 


bound<>.ry  where 

=  E^°^(k-K)  (46) 

we  examine  the  pole  structure  a  little  more  carefxilly.  One  then  neglects  all 
other  terms  ercept  this  one  in  (43)  so  that  the  pole  is  now  determined  by 
the  equation 


(w-E^°^k))  (w-E^^^k-K))  - 


V^^V) 

ciKr" 


=  0 


(47) 


If  we  ignore  the  <j  dependence  in  M  ,  we  arrive  at  a  familiar  looking 


u  = 


E^^^(k)4E^^^(k-K)  I  ^ 


edge 

//E^°^(k)  -E^®^(k-K)^ 

1  ^  1  (l)  1  ^ 

A  ^  j 

^  1  €(K)  1 

'\  1 

(48) 


This  generalizes  the  nearly  free  electron  gap  equation  at  the  zone  edge  [6] 
for  the  interacting  electron  system.  Recently,  Tanaka  [7]  has  also  examined 
this  problem  using  diagrammatic  approach. 


DISCUSSION 

We  have  here  constructed  a  perturbation  method  for  self-consistent 
equations  in  a  Green's  function  formulation.  This  is  a  generalization  of 
Peng's  work  for  static  HE  equations.  We  incidentally  show  that  the  amplitude 
of  the  first-order  wave  function  satisfy  the  usual  RPA  equations.  As  an 
example  of  the  method  we  briefly  outline  a  calculation  of  the  effect  of 
periodic  field  on  a  many-electron  system. 
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